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ABSTRACT 

Preferential solvation analysis on experimental solubility data 
of moxidectin was performed by using the Inverse Kirkwood- 
Buff Integrals approach at 298.15 K. Local mole fraction of 
solvents in the solvation shell of moxidectin as well as the 
extent of preferential solvation by each of solvents were 
calculated as a function of the bulk mole fraction of binary 
mixtures of water with methanol, ethanol, isopropanol, and 
ethylene glycol. Results indicate that preferential solvation of 
moxidectin by water occurs in water-rich regions. Whereas, 
beyond water-rich regions, moxidectin was preferentially 
surrounded by alcohols in all binary mixtures studied. 
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1. Introduction 

Investigation of intermolecular interactions in multicomponent systems and their relationship with the system 
structure has always been of interest to researchers [1-5]. As an attractive issue in solution chemistry, preferential 
solvation may occur in mixtures due to differences in interactions between components. In such cases, components 
experience different compositions in their solvation shell with respect to the bulk. The Inverse Kirkwood-Buff 
Integrals (IKBI) provides the best applicable approach for the analysis of preferential solvation of solutes in a mixture 
of two solvents [6]. In this method, the probability of finding a molecule of one solvent around a solute molecule is 
defined by the integrals of radial distribution functions for the pair of two molecules. With the help of statistical 
thermodynamics in an inverse way, KBI can be determined from some macroscopic thermodynamics data for mixed 
solvents and solutions. Analysis of obtained KBI over the mixture composition gives information on microscopic local 
composition of solvation shell of components and thus on the probability and extent of preferential solvation [1-6]. 

Solvent screening is the main step in process of optimization of synthesis, crystallization, extraction, and 
formulation of drugs in the pharmaceutical industry. In this regard, knowledge of solubility and dissolution 
mechanisms is critical. Most drugs because of their organic nature show low aqueous solubility that restricts their 
pharmacokinetics and pharmacodynamics. Co-solvency method in which miscible organic solvents are mixed with 
water is an advisable convenient and effective approach for solubility enhancement [7]. Despite of widely used co- 
solvency method in pharmacy, investigation on preferential solvation aimed to understand how a change in mixture 
composition improves the solubility has only been carried out in recent years [4, 8-11]. 

Moxidectin (with the molecular structure shown in Figure 1) is used as a drug in the treatment of various diseases 
in a veterinary clinics. Recently, the solubility of moxidectin has been determined in aqueous mixtures of alcohols 
including methanol, ethanol, isopropanol, and ethylene glycol over a temperature range of 278.15—323.15 K [12]. 
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Figure 1. Molecular structure of moxidectin 


Different co-solvency models have been used to correlate data to being used for the prediction of solubility over 
other compositions. However, no further analysis has been done on experimental data in the binary mixture in order 
to gain insight into solvation and enhance the effect of co-solvent on solubility. Most studies on drug solubility data 
in mixed solvents have been devoted to mathematical modeling and possibly predicting solubility from drug solubility 
in pure solvents. Depending kind of solvent mixtures, models such as Jouyban-Acree, NRTL, Wilson, and UNIQUAC 
have been developed [13-16]. Such models are successful to make a correlation between solubility data with mixture 
composition and temperature, but provide little detail on the actual molecular environment of the solute and the 
interactions leading to preferences between solute and each solvent in the mixture. The Inverse Kirkwood-Buff 
Integrals (IKBI) model is a powerful approach for providing insight into microscopic structures in mixed components 
based on thermodynamic data analysis [17-23]. Therefore, the aim of this work is preferential solvation analysis of 
solubility data of moxidectin in a binary mixture with the help of the IKBI approach to determine the local mole 
fractions of solvent in the solvation shell of moxidectin as well as the extent of preferential solvation over the mixture 
composition. 


2. Theoretical considerations 
2.1. IKBI model 
The quantity Xi s= xis — Xi explains the preferential solvation of solute s by solvent i. Here, xis is the local mole 
fraction of solvent i around the solute s, and x; is the mole fraction of solvent i in the bulk mixture of two solvents. 
According to Kirkwood-Buff Integrals, the extent of preferential solvation of can is calculated by Eq. (1) [17-23]. 
XXa (Gs =Gy;) (1) 
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In this equation, an affinity between components is defined by Gi s so-called Kirkwood-Buff integral. From point of 
view of statistical thermodynamics, Gis can be readily determined by the following integration on the radial 
distribution function between solute and solvent, gi,s(r) as Eq. (2). 


_ Toor 2 2 
G, = f (g,,-D4ar°dr (2) 
here, r is the distance between centers of moiccules. However, Gis can be accurately determined in an inverse 
approach from the macroscopic properties of the solution by using Eq. (3) and (4). 


G,, =RTK,-V, Ped 
Q (3) 


G, =RTk, -V, anp 
£ (4) 

In these equations, isothermal compressibility xr molar volume for organic solvent (methanol, ethanol, isopropanol, 
and ethylene glycol) Vi, molar volume for water V2, and molar volume for moxidectin V, are needed to be known. R 
and T are the universal gas constant and temperature, respectively. Other functions D and Q are determined by 


thermodynamic functions including excess Gibbs free energy of mixing of two solvents in a binary mixture, Ga and 
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the standard molar Gibbs energies of solute’s transfer, AGs 22 from pure solvent 2 to binary mixture by using the 
following equations Eq. (5) and (6). 


Gr 
O=Rr ss 5 J (5) 
%2 Ir p 
9A, Ge 
D = tr ~ (s,2314+2) 
l dx, l P (6) 


Binary mixtures of two solvents AGs 2n) can be readily determined by the mole fraction solubility of solute xs, 
based on Eq. (7). 
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Finally, in Eq. (1), the correlation volume Voor is the area around the solute in which preferential solvation occurs. 
It is obvious that Vior depends on the local composition. Eq. (10) is recommended for the calculation of Veo. by 


considering the first shell around the solute [4]. 
3 


a 2522.5) +0.1363( 2. +x4.V,)" -0.085] (8) 


Therefore, one can analyze the preferential solvation by solving Eq. (1) by prior simultaneous solving of Eqs. (3)- 
(8) by knowledge of thermodynamic properties of solution [17-25]. 


3. Results and discussion 

For calculation of function Q in Eq. (5), the excess Gibbs free energy of mixing water with methanol, ethanol, 
isopropanol, and ethylene glycol was extracted from the reference at 298.15 K [26]. Partial molar volume for solvents 
in mixtures was calculated from density data at 298.15 K available for the studied binary mixture from the literature 


[27, 28]. The isothermal compressibility for each aqueous binary mixture can be safely computed by x1 ATI+ 


X22. for pure solvents, corresponding data was taken from the literature [29]. The standard Gibbs free energy of 
transfer of moxidectin from water to methanol (1) + water (2), ethanol (1) + water (2), isopropanol (1) + water (2), 
and ethylene glycol (1) + water (2) mixtures was determined by Eq. (7). Data for the mole fraction solubility of 
moxidectin in binary mixtures were taken from the literature [12]. In addition functions, D, Q, Gis, and G2, were 
calculated by Eqs. (3)-(6) versus the bulk mole fractions and tabulated in Table 1-4. The correlation volume in Eq. 
(8) is a function of the local mole fractions. Therefore, in an iterative process, dx1,s, and Veor were computed by 
simultaneous solving of Eqs. (1) and (8) through iteration on the local mole fractions until a non-variant value was 
obtained. 
Table 1. Kirkwood-Buff functions D, Q, G2,s and Gi,s in methanol (1) + water (2) mixture at 298.15 K. 


D Q Gis G25 
(kJ-mol=') (kJ-mol=') (cm3.mol!) (cm?.mol'!) 
0.00 -194.04 2.48 -1929 -516 
0.05 -97.71 2.51 -1184 -589 
0.10 -58.68 2.47 -902 -604 
0.15 -39.10 2.38 -768 -608 
0.20 -27.91 2.27 -693 -608 
0.25 -20.91 2.14 -647 -608 
0.30 -16.25 2.01 -616 -608 
0.35 -12.99 1.90 -594 -608 
0.40 -10.62 1.81 -577 -607 
0.45 -8.85 1.74 -563 -605 
0.50 -7.48 1.71 -552 -601 
0.55 -6.41 1.72 -543 -597 
0.60 -5.55 1.76 -536 -591 
0.65 -4.86 1.83 -530 -584 
0.70 -4.29 1.93 -525 -577 
0.75 -3.81 2.04 -522 -571 
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0.80 -3.41 2.17 -519 -565 
0.85 -3.07 2.29 -517 -560 
0.90 -2.711 2.39 -516 -556 
0.95 -2.52 2.46 -515 -554 
1.00 -2.30 2.48 -514 -552 


The affinity between alcohol and moxidectin, and between water and moxidectin are respectively determined by 
Kirkwood-Buff integrals G1, and G72,. It is evident from Tables 1-4 that G1 s, and G2, were negative in all compositions 
for all aqueous binary mixtures. It means that moxidectin show affinity for positive interaction with both solvents in 
all binary mixtures studied. However, a comparison of their values indicates that in water-rich regions in all binary 
mixtures, G2, is higher than Gs. This result indicates that moxidectin has a more affinity for solvation by water with 
respect to alcohols in water-rich regions. 


Table 2. Kirkwood-Buff functions D, Q, G2,s and Gi,s in ethanol (1) + water (2) mixture at 298.15 K. 


A D Q Gis G25 
(kJ-mol"') (kJ-mol"') (cm3.mol!) (cm?.mol!'!) 
0.00 -336.84 2.48 -2969 -516 
0.05 -118.90 2.22 -1434 -660 
0.10 -60.13 2.04 -994 -676 
0.15 -36.20 1.91 -805 -672 
0.20 -24.15 1.81 -705 -664 
0.25 -17.25 1.73 -647 -655 
0.30 -12.94 1.66 -611 -648 
0.35 -10.06 1.58 -587 -643 
0.40 -8.05 1.49 -570 -640 
0.45 -6.58 1.39 -559 -639 
0.50 -5.48 1.27 -551 -640 
0.55 -4.64 1.15 -545 -644 
0.60 -3.97 1.04 -540 -649 
0.65 -3.44 0.94 -536 -654 
0.70 -3.01 0.87 -531 -656 
0.75 -2.66 0.86 -527 -650 
0.80 -2.36 0.92 -523 -635 
0.85 -2.11 1.08 -519 -612 
0.90 -1.90 1.37 -517 -588 
0.95 -1.72 1.82 -515 -567 
1.00 -1.56 2.48 -514 -551 


Table 3. Kirkwood-Buff functions D, Q, G2,s, and Gi,s in isopropanol (1) + water (2) mixture at 298.15 K. 


s D Q Gis G25 
(kJ-mol"!) (kJ-mol"!) (cm?.mol!) (cm?.mol t) 
0.00 -230.71 2.48 -2196 -516 
0.05 -94.26 2.40 -1187 -654 
0.10 -50.91 2.26 -878 -677 
0.15 -31.80 2.07 -746 -684 
0.20 -21.74 1.85 -680 -689 
0.25 -15.79 1.61 -643 -699 
0.30 -11.99 1.36 -621 -715 
0.35 -9.41 1.11 -608 -739 
0.40 -7.58 0.88 -602 -776 
0.45 -6.24 0.68 -600 -831 
0.50 -5.23 0.51 -600 -906 
0.55 -4.44 0.38 -601 -998 
0.60 -3.82 0.32 -595 -1069 
0.65 -3.32 0.31 -577 -1050 
0.70 -2.91 0.37 -554 -938 
0.75 -2.57 0.50 -535 -807 
0.80 -2.29 0.72 -525 -709 
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0.85 -2.06 1.02 -519 -645 
0.90 -1.85 1.41 -516 -604 
0.95 -1.68 1.90 -515 -578 
1.00 -1.53 2.48 -514 -561 


Table 4. Kirkwood-Buff functions D, Q, G2,s, and Gi,s in ethylene glycol (1) + water (2) mixture at 298.15 K. 


D Q Gis G25 
(kJ-mol"') (kJ-mol"') (cm3.mol!) (cm?.mol'!) 
0.00 -169.10 2.48 -1747 -516 
0.05 -81.44 2.56 -1062 -602 
0.10 -47.28 2.61 -811 -615 
0.15 -30.69 2.63 -695 -612 
0.20 -21.46 2.63 -633 -606 
0.25 -15.81 2.63 -597 -599 
0.30 -12.12 2.63 -574 -593 
0.35 -9.58 2.63 -558 -587 
0.40 -7.76 2.63 -548 -582 
0.45 -6.41 2.64 -540 -577 
0.50 -5.38 2.66 -534 -572 
0.55 -4.58 2.69 -530 -568 
0.60 -3.95 2.72 -526 -565 
0.65 -3.44 2.76 -524 -561 
0.70 -3.02 2.79 -522 -559 
0.75 -2.67 2.82 -520 -556 
0.80 -2.38 2.82 -519 -554 
0.85 -2.13 2.80 -518 -552 
0.90 -1.93 2.75 -517 -552 
0.95 -1.75 2.64 -517 -551 
1.00 -1.59 2.48 -516 -552 


The extent of preferential solvation of moxidectin by alcohol OX1,; was determined as a function of mole fraction 
of alcohol in aqueous binary mixtures, and plotted against mole fraction of organic solvent in Figure 2. 6x1,; was 
negative at 0.00 < x; < 0.15 for isopropanol (1) + water (2), at 0.00 < xı < 0.20 for ethanol (1) + water (2), and ethylene 
glycol (1) + water, at 0.00 < xı < 0.30 for methanol (1) + water (2) mixtures. Therefore, moxidectin was preferentially 
solvated by water in water-rich regions in all binary aqueous mixtures studied here. The reason for this observation 
may be because of the enhancement effect of alcohols on the structuring of the hydrogen bonding water network 
around the solute [30]. 


Figure 2. Preferential solvation of moxidectin by organic solvent as a mole fraction of methanol (©), ethanol (0), isopropanol (e), 
and ethylene glycol (A). 
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In water-rich region, the order of preferential solvation by water in mixtures is as ethanol (1) + water (2) > methanol 
(1) + water (2) > isopropanol (1) + water (2) > ethylene glycol (1) + water (2). Marcus showed that water molecules 
become more structured in presence of organic solvents such as alcohols. However, the enhancement effect was the 
least for ethylene glycol [31]. Furthermore, moxidectin can interact with water as a Lewis base by its oxygen atoms. 
Water is more acidic than alcohols and thus is more favorable to making such hydrogen bonding interactions with 
moxidectin. The 6x); was positive for all binary mixtures beyond water-rich regions up to pure alcohol. This result 
reveals that moxidectin was solvated more selective by methanol, ethanol, isopropanol, and ethylene glycol rather 
than by water in these regions. Maximum value for preferential solvation was observed at x;=0.55 for methanol (1) + 
water (2), at x1=0.60 for ethanol (1) + water (2) and isopropanol (1) + water (2), and at x,;=0.45 for ethylene glycol (1) 
+ water (2) mixtures. The reason for preferential solvation by alcohols may be related to the strong network structure 
of water. Cavity formation in water is a much more energy-consuming step in the solvation of moxidectin with respect 
to alcohols. In addition, all alcohols methanol, ethanol, isopropanol, and ethylene glycol are more basic than water, 
and thus can make more favorable interactions with acidic hydroxyl sites on moxidectin. It seems that moxidectin 
prefers to make such hydrogen bonding interactions with alcohols that with water in aqueous binary mixtures studied. 
The order of preferential solvation by alcohol, 6x15, in these regions is as isopropanol (1) + water (2) > ethanol (1) + 
water (2) > methanol (1) + water (2) > ethylene glycol (1) + water (2). This observation is in good agreement with the 
order of hydrogen bond acceptor capacity of solvents described by the Kamlet-Taft hydrogen bond acceptor parameter 
as $= 0.84 for isopropanol > f= 0.75 for ethanol > = 0.66 for methanol > 2 = 0.52 for ethylene glycol > B= 0.47 for 
water [32]. It confirms that the hydrogen bonding interactions in which moxidectin acts as Lewis acid are mainly 
responsible for preferential solvation in these regions. 


4. Conclusion 

In this work, the approach of Inverse Kirkwood-Buff Integrals was successfully used for the analysis of 
experimental data of moxidectin in binary mixtures of water with methanol, ethanol, isopropanol, and ethylene glycol. 
Gibbs free energy of transfer of moxidectin from water to the mixture was calculated and used with other 
thermodynamic properties of mixture for computation of local mole fraction of solvents around moxidectin and extent 
of preferential solvation. Kirkwood-Buff integral functions and preferential solvation parameter indicate that 
moxidectin was solvated more by water in the water-rich region, and by alcohols in other remaining compositions up 
to pure organic solvent. 
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